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It was long believed that the singularity inside a realistic, rotating black hole must be spacelike. 
However, studies of the internal geometry of black holes indicate a more complicated structure is 
typical. While it seems likely that an observer falling into a black hole with the collapsing star 
encounters a crushing spacelike singularity, an observer falling in at late times generally reaches 
a null singularity which is vastly different in character to the standard Belinsky, Khalatnikov and 
Lifschitz (BKL) spacelike singularity [V. A. Belinsky, I. M. Khalatnikov, and E. M. Lifshitz, Sov. 
Phys. JETP 32, 169 (1970)]. In the spirit of the classic work of BKL we present an asymptotic 
analysis of the null singularity inside a realistic black hole. Motivated by current understanding of 
spherical models, we argue that the Einstein equations reduce to a simple form in the neighborhood 
of the null singularity. The main results arising from this approach are demonstrated using an almost 
plane symmetric model. The analysis shows that the null singularity results from the blueshift of 
the late-time gravitational wave tail; the amplitude of these gravitational waves is taken to decay 
as an inverse power of advanced time as suggested by perturbation theory. The divergence of the 
Weyl curvature at the null singularity is dominated by the propagating modes of the gravitational 
field, that is 

ri rvxp-tS , T , , T , ~{21+3)2kv 

as v — » oo at the Cauchy horizon. Here, and ^4 are the Newman-Penrose Weyl scalars, and 
/ > 2 is the multipole order of the perturbations crossing the event horizon. The null singularity 
is weak in the sense that tidal distortion remains bounded along timelike geodesies crossing the 
Cauchy horizon. These results are in agreement with previous analyses of black hole interiors. We 
briefly discuss some outstanding problems which must be resolved before the picture of the generic 
black hole interior is complete. 

Pacs numbers: 03.70.+k, 98.80.Cq 
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I. INTRODUCTION 

Spacetime singularities arc an inevitable consequence 
of the Einstein field equations. They mark the boundary 
of spacetime, and the limit of our current understand- 
ing of gravitational physics. Unfortunately, the powerful 
techniques used to demonstrate the existence of singu- 
larities say nothing about the character of these physical 
blemishes 0. 

The classic works of Belinsky, Khalatnikov and Lifs- 
chitz || address this deficiency by integrating Einstein's 
equations in the neighborhood of a spacelike singular- 
ity. They present compelling evidence that the general 
solution takes an inhomogeneous Kasner form exhibit- 
ing chaotic oscillations of the Kasner axes as a crushing 
singularity is approached 0. The functional genericity 
of the BKL solutions is widely interpreted as indicating 
that all physical singularities must be of this form. 

Nonetheless, studies of the internal geometry of black 
holes suggest a picture involving two distinct regimes. 
Observers falling into the black hole with the collapsing 
star generally encounter a spacelike singularity (which 
is presumably of the BKL type). On the other hand, 
observers that fall in at late times, when the external 
geometry has settled down to an almost stationary state, 
encounter a weak, null singularity]^] of a type similar to 
the mass- inflation singularity of Poisson and Israel || . 

To understand this behavior we must first discuss the 
formation of a black hole by the gravitational collapse of 
a rotating star. At late times, the external gravitational 
field is believed to settle down to a Kerr-Newman solu- 
tion. These solutions have a timelike singularity which 
is preceded by a Cauchy horizon — a null hypersurface 
marking the boundary of the domain of dependence for 
Cauchy data prescribed in the black hole exterior. The 
Cauchy horizon is non-singular, and the spacetime can 
be analytically extended through it. The global solu- 
tion then suggests that black holes act as tunnels from 
our asymptotically flat universe to other identical, but 
distinct, universes. Indeed there is an infinite lattice of 
universes extending into the past and future of our own. 
Furthermore, observers may travel through the ring sin- 
gularity inside these black holes, passing to achronal re- 
gions of spacetime. 

In the late 1960's Penrose pointed out that the Cauchy 
horizon inside such a black hole is unstable |6). Time- 
dependent perturbations originating outside the black 
hole get infinitely blucshiftcd as they propagate inwards 
near to the Cauchy horizon, consequently, the energy 
density associated with these perturbations diverges as 
measured by a free-falling observer attempting to cross 
through the horizon. Perturbative calculations M in 



*The weakness of the mass-inflation singularity was first elu- 
cidated by Ori in Refs. 



Reissner-Nordstrom and Kerr spacetimes have validated 
Penrose's original arguments. 

In general, gravitational collapse is expected to be 
asymmetric, so that gravitational and electro-magnetic 
waves arc emitted by a newly formed black hole as it set- 
tles down to a stationary, axisymmetric state. Detailed 
studies of perturbations of black hole geometries show 
that such gravitational wave emission results in wave tails 
which decay according to an inverse power law of time in 
the exterior of the black hole Some of this radia- 

tion inevitably crosses the event horizon getting infinitely 
blueshifted near to the Cauchy horizon. For this reason 
one might expect the internal geometry of a black hole 
formed by collapse is significantly different to that of the 
exact stationary solutions. 



A. Spherical models of black hole interiors 

The first attempts to understand the back-reaction 
of the blueshifted, radiative tail on the internal geom- 
etry of black holes were restricted to spherical symmetry. 
Hiscock |lCj argued that Isaacson's |lf|] effective stress- 
energy description of high-frequency gravitational waves 
should be valid near the Cauchy horizon. He considered a 
charged black hole with a directed influx of lightlike dust 
with stress-energy tensor T a p = Pi n (. a £p where i a l a = 0, 
and showed that an observer dependent singularity forms 
along the Cauchy horizon in this circumstance. 

In reality, some of the infalling radiation is back- 
scattered off the curvature inside the event horizon. Pois- 
son and Israel j3| modeled this effect by another flux of 
lightlike dust moving to the right; they demonstrated, 
for the first time, that non-linear effects transform the 
Cauchy horizon into a scalar curvature singularity. It is 
worth summarizing the essence of their argument here. 

Figure [l] shows the setting for the characteristic initial 
value problem. Charged ingoing and outgoing Vaidya 
solutions in regions II and III, respectively, are matched 
continuously onto region I, which is described by a static 
Reissner-Nordstrom solution with mass too and charge q. 
(The pure ingoing region, II, was studied by Hiscock.) In 
region IV the line element can be written as 



ds 2 



-2e x dudv + r 2 (d6 2 



(1.1) 



where r = r(u,v) and A = \{u,v). The coordinates are 
chosen so that v is standard advanced time, i.e. v = oo 
at future null infinity, and the retarded time u = — oo 
on the black hole event horizon. The stress-energy of 
cross-flowing null dust is 
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t(«). 
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(1.2) 



where i a = —d a v and n a = —d a u. The luminosity 
function L- ln (v) is fixed by requiring that the flux of 
stress-energy across the event horizon decays as an in- 
verse power of advanced time, thus 
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L in (v) = a(Kv)~ 



(1.3) 



with q = Al + 6, where / is the multipole order of the per- 
turbing field, fixed by Price's analysis The constant 
a depends on the luminosity of the star that collapses 
to form the black hole, and k is taken to be the surface 
gravity of the stationary segment of the Cauchy horizon 
in region II. This functional form is motivated by our un- 
derstanding of radiative tails in the exterior of the black 
hole |f| . The outflux is produced by scattering of ingoing 
tail radiation inside the black hole, and consequently also 
has a power law form |7p^-p^| 



Lout(u) = I3(~ku) q 
for large negative values of the coordinate u. 



(1.4) 




FIG. 1. The Poisson-Israel spherical model of the black 
hole interior. The infalling wave-tail of gravitational radia- 
tion, and its scattered component, are modeled by lightlike 
dust. The region of purely ingoing (outgoing) dust is indi- 
cated by II (III) in the figure, and is described by a charged 
Vaidya solution. The solutions in regions II and III are 
matched to a Reissner-Nordstrom solution (region I) along 
the null rays v = vq and u — uo . As usual I + and T~ denote 
future and past null infinity. The inflow causes the apparent 
horizon (AH) to expand and smoothly approach the event 
horizon (EH) of the final black hole. The outflow in region IV 
causes the, initially static, Cauchy horizon (CH) to contract. 
This, together with the fatal blue shift experienced by the 
inflow, causes a curvature singularity to develop along CH. 

In the limit as v — > oo, the corresponding solution of 
Einstein's equations is well approximated by 



' CH 



+ 



2a{nv)- q+1 
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K 2 (p — 1) 



(1.5) 
(1.6) 



where tch is the constant radius of the Cauchy horizon in 
region II. The square of the Weyl tensor diverges on the 
Cauchy horizo n in this solution, yet the radial function 
given by Eq. (1.6) is non-zero for sufficiently large \u\. 
This indicates that the singularity is not a central (r = 0) 
singularity with which we are familiar, in fact it is a null 
singularity. Spherical symmetry allows the introduction 
of a mass function m(u, v) which is directly related to the 
Weyl curvature of the spacetime: 



m(u,v)= (iC7 e %r 3 +g 2 A) 
afie KV 



(Kv)- q (-Ku)- q 



(1.7) 
(1.8) 



The divergence of this mass function as v — > oo in region 
IV prompted Poisson and Israel to refer to the accom- 
panying scalar-curvature singularity as a mass inflation 
singularity. 

Subsequently, Ori Q pointed out that physical objects 
which encounter a mass-inflation singularity experience 
only finite tidal distortion, even though the tidal forces 
acting on the object diverge. In this sense, the mass- 
inflation singularities are weak alth oug h the relevance of 
this point continues to be debated Jlqjl6| . 

A variety of spherical models have now been stud- 
ied flj7|-fl9|| , all of them indicate the presence of a null, 
scalar curvature singularity at the Cauchy horizon. The 
numerical analyses in |l^] also demonstrates the slow con- 
traction of the null generators of the Cauchy horizon to 
zero radius, and the formation of a spacelike singularity 
deep inside the black hole core. 



B. Beyond spherical symmetry 

While simplified models argue strongly in favor of 
weak, null singularities inside black holes, the results are 
restricted entirely to spherical symmetry. Nevertheless 
the physics behind the mass-inflation phenomenon is ex- 
tremely general. Perturbations originating in the exter- 
nal universe get infinitely blueshifted as they propagate 
close to the Cauchy horizon resulting in a scalar curva- 
ture singularity through non-linear interaction with grav- 
ity. Attempts have been made to investigate the problem 
in less symmetric situations. Ori j^] has used non-linear 
perturbation theory to examine back-reaction of pertur- 
bations on the Kerr geometry. His results support the 
existence of a weak, null singularity in this case. Bo- 
nanno (20j matched two radiating Kerr solutions along a 
null hypersurface and showed, in the limit of small angu- 
lar momentum, that the mass-function diverges on a null 
hypersurface. 

Brady and Chambers demonstrated that weak, 
null singularities are consistent with the Einstein equa- 
tions on a pair of intersecting null surfaces, and have 
investigated the rate of divergence of curvature on a 
null surface crossing the singularity. An important step 
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has also been taken by Ori and Flanagan [g2| who ar- 
gue that the vacuum, Einstein equations admit function- 
ally generic solutions containing weak, null singularities. 
Their results invalidate local arguments Q that null 
singularities are always unstable to transformation into 
crushing spacelike singularities of the BKL type. 



C. Discussion and overview of this paper 

In this paper we demonstrate that a null singularity 
replaces the Cauchy horizon inside a black hole formed 
by gravitational collapse. The method employed in the 
following analysis is similar in spirit to that of BKL [^) ; 
we use an asymptotic expansion to find the leading order 
behavior of the metric near the Cauchy horizon with- 
out making any assumptions about the symmetry of the 
spacetime. In this way we show that the wave tail of 
gravitational collapse results in the formation of a null 
curvature singularity at the Cauchy horizon. 

Our results should be compared with those obtained 
using non-linear perturbation theory. In that context, 
Ori H has previously argued that a weak, null singular- 
ity is present along the Cauchy horizon of a generic, ro- 
tating black hole. Our analysis, being asymptotic, relies 
on physical arguments to provide initial data near to the 
Cauchy horizon; in this sense, it is a local analysis. Nev- 
ertheless, the results presented below embody all orders 
of perturbation theory and lend further support to the 
claim that the perturbative approach captures the essen- 
tial features of spacetime structure near the Cauchy hori- 
zon singularity inside a realistic, rotating black hole 

It is natural to use null coordinates to attack this prob- 
lem since the Cauchy horizon is a null hypersurface. To 
facilitate the analysis we employ a double-null formalism 
pij in which spacetime is decomposed into two families 
of intersecting null hypersurfaces. This decomposition of 
the Einstein field equations is reviewed in Sec. |ll] where 
it is shown that the most general double-null metric de- 
pends on six functions of four variables. Interested read- 
ers are referred to Brady et al for more details. 

The main assumption made in this paper is that the 
results of scattering on the Kerr background |2[| are suf- 
ficient to determine the initial data for the interior prob- 
lem. This assumption and other physical considerations 
are discussed in Sec. 



Ill 



In Sec. IV we present the essential features of our argu- 
ments in a simplified context where the gravitational field 
has only one dynamical degree of freedom. We carefully 
construct the solution in the neighborhood of the Cauchy 
horizon for this model problem. Moreover, we demon- 
strate that the solution is characterized by a singularity 
at which the Weyl curvature diverges like 



(1.9) 



the singularity, where it becomes asymptotically Petrov 
type N. 

The relevance of the simplified model is made clear 
in Sec. where we present the general analysis. By 
neglecting terms in the field equations which we expect 
to be exponentially small near to the Cauchy horizon, it 
is shown that the general solution is identical in character 
to that presented in Sec. |y|. 

Finally, we discuss the strength of the singularity in 
Sec. fv|. While tidal forces experienced by an observer 
diverge at the singularity, the integrated tidal distortion 
is finite. As discussed by Ori the singularity is 

therefore weak and we are once again faced with the pos- 
sibility that spacetime might extend beyond the Cauchy 
horizon (in a non-analytic manner). We emphasize that 
gravitational shock waves are the only type of curvature 
which can be confined to a thin layer in the absence of 
matter This suggests that a classical continua- 

tion of spacetime beyond the Cauchy horizon is unlikely 

We adopt Misner, Thorne and Wheeler curvature con- 
ventions p7[ with metric signature ( — V H — V) through- 
out the text. Greek indices a, f3, . . . run from to 3; 
upper-case Latin indices A, B, . . . take values (0, 1); and 
lower-case Latin indices a,b, . . . take values (2,3). 



II. DOUBLE NULL FORMALISM 

The complexity of the Einstein equations requires the 
careful selection of a formalism with which to pursue our 
goal of understanding the nature of spacetime near null 
singularities. In |24|] we have developed the necessary ma- 
chinery based on a dual null decomposition of the Ein- 
stein field equations. In this section we present those 
details which are necessary for the application at hand. 
The interested reader is referred to |24|] for a more com- 
plete treatment. 



A. Lightlike foliation of spacetime 

We suppose that we are given a foliation of spacetime 
by lightlike hypersurfaces E° with normal generators £^ , 
and a second, independent foliation by lightlike hyper- 
surfaces E 1 with generators £W nowhere parallel to £^°\ 
The intersections of {E } and {E 1 } define a foliation of 
codimension 2 by spacelike 2-surfaces S. (The topology 
of S is unspecified. All our considerations are local.) S 
has exactly two lightlike normals at each of its points, 
co-directed with £^ and £^\ 

In terms of local charts, the foliation is described by 
the embedding relations 



x a =x a (u A ,9 a ). 



(2.1) 



Here. 



as v — > oo at the Cauchy horizon. We show that the 
resulting spacetime is not algebraically special except at 



are four-dimensional spacetime co-ordinates; 



u = v and u 1 = u are a pair of scalar fields constant over 
each of the hypersurfaces E° and E 1 respectively. The 
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intrinsic co-ordinates 8 a of the 2-spaces S 1 , each charac- 
terized by a fixed pair of values (u ,^ 1 ), are convected 
along the vectors d/du A . 

It is convenient to effect a partial normalization of the 
lightlike vectors by imposing the condition 

WgoPfW^B) = e x j]AB _ e A anti _ diag (_ 1; _ 1 ) (2.2) 

for some scalar field X(x a ), reflecting the freedom to ar- 
bitrarily rescale a null vector. The superscript (4) is used 
to indicate four dimensional objects whenever confusion 
might arise. We can use r\ AB and its inverse t\ab to raise 
and lower uppercase Latin indices. Furthermore, the vec- 
tors 6- A ^ are parallel to the gradients of u A allowing us 
to write 

6 A ) = e x d a u A . (2.3) 

The pair of vectors e"^ , defined from ( fO| ) by 

e a (a) = dx a /d9 a , (2.4) 

are holonomic basis vectors tangent to 5*. The intrinsic 
metric g ab of S is determined by their scalar products: 



(4) a 

9ab = { ' ' g a pe {a) e p {b) . 



(2.5) 



Lower-case Latin indices are lowered and raised with g ab 
and its inverse g ab . Since £( A ' is normal to every vector 
in T, A , we have 



£ (A) e? n = 



(2.6) 



Finally, introducing a pair of shift vectors tangent 
to S by 



dx a 



(2.7) 



an arbitrary displacement dx a in spacetime is 

dx a = l a [A) du A + ef a) (d9 a + s a A du A ) . (2.8) 

Thus, the final form of the line element used below to 
discuss the black hole interior is 

ds 2 = e x r lAB du A du B + g ab {dd a + s a A du A ){dO b + s b B du B ). 

(2.9) 

We will impose further coordinate conditions as the need 
arises, but it is manifest that Eq. (2.9) is sufficient to de- 
scribe the most general spacetime admitted as a solution 
to the Einstein field equations. 



must first introduce some notation. The extrinsic curva- 
tures of S are defined by 



2K Aab — e™ a )e^ 



(4)£ 



l (A) 9af3 



DA9ab , 



(2.10) 



£a ep g a b is tne 



where g a p = {4, g a - e x ri AB t a 1^ 
spatial projection tensor, and ^Ce (A) is the four dimen 
sional Lie derivative along l/j^ . We have also introduced 
the derivative operators D A acting on two-dimensional 
tensorial objects according to the rule 



D A X a - 



b... 



(2.11) 



Here C«d is the two-dimensional Lie derivative. (Gener- 

A . . 

ally D a can be defined in a four-dimcnsionally covariant 

manner by its action on spatial four tensors pjfl; 



D A X° 



,(a) P P (4), 



r^>e lb) ^C kA) {X c --- d ..e^), (2.12) 



thus providing a complete geometrical interpretation of 
Da in terms of Lie derivatives.) 

The geometric meaning of the extrinsic curvature is 
illustrated by Lie propagating a circle on S along £(a)- 
The expansion rate of the light rays is given by the trace 
of the extrinsic curvature Ka = g ab KAab- The traceless 
part of the extrinsic curvature, UAab = KAab — ^gabKA, 
is the shear of the null congruence, and the twist is 



af \VBS A - S B d b S A ) 



(2.13) 



where e AB is the anti-symmetric tensor density (eoi = 1)- 
Clearly uj a is purely spatial. 

Finally, the expressions for the Ricci tensor are 



B. 2+2 decomposition of curvature 

Following |24|] , the components of the Ricci tensor can 
be compactly expressed in terms of two dimensionally 
covariant quantities. In passing to these expression we 
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^Rab — R, 



a (3 
*0 e (a) e (b) 



Rab — RasP^Z 



(Af(B) 



RAa - Raf3ef a) l^ A) 



-y ab e- x (D A K A + K A K A ) + e- x (-D A a Aab + 2a A(a d a A b)d ) 

2 (2)R 9ab - \ e ' 2X ^a^b - \ab - T^A a A j6 

-D/ A K B ) — \K A K B — o Aab a B a + K( A D B )\ 



—$\ab [{D E + K E ) D E X - e- X Lu a u; a + (e A ) ;Q ;a ] 
= cTAa b - b - \d a K A - ^d a D A X + ^K A d a X 

+ \e AE e- X [{D B + K B ) uj a - uj a D B \] , 



(2.14) 
(2.15) 

(2.16) 



where ( 2 >R is the curvature scalar associated with 
the two- metric g ab , and a semi-colon indicates the two- 
dimensional covariant derivative. 



£ and near to the Cauchy horizon, given initial data con- 
sistent with scattering of a test gravitational wave field 
inside a Kerr black hole [Eal. 



III. ASSUMPTIONS 

The mathematical theory of black holes is well estab- 
lished |28| . Indeed, it is widely accepted that the external 
geometry of a black hole formed by the gravitational col- 
lapse of a star is described by a Kerr-Newman solution 
at late times. The mechanism by which the black hole 
settles down to this stationary state was elucidated by 
Price [|| for nearly spherical, uncharged collapse, and his 
results have been extended to other situations by several 
authors ]g|j?9[]. These works provide the initial condi- 
tions for the interior problem by determining the behav- 
ior of the radiative tail of gravitational waves at the event 
horizon of the black hole: in particular, the amplitude of 
the gravitational wave flux generally decays as an inverse 
power of advanced time. 

Some simplifying assumptions are required in order to 
make progress against the black-hole interior problem. 
We clarify our approach below. 



A. Physical considerations 

The inner Cauchy horizon of a stationary black hole 
is unstable to linear perturbations that undergo infinite 
gravitational blueshift there. The central assumption of 
our analysis is that an asymptotic limit exists in which 
effects on the geometry, inside a black hole formed by 
gravitational collapse, are dominated by gravitationally 
blueshifted tail-radiation which propagates into the black 
hole at late times (as measured by external observers). 
The postulated structure of spacetime inside the black 
hole is shown in Fig. E|. On physical grounds |3(| one 
expects that the region between the event horizon and 
the three dimensional spacelike surface E can be ade- 
quately described by linear perturbation theory — there 
is no physical mechanism to excite an instability of the 
interior before the blueshifting takes hold. We examine 
the structure of spacetime in the region to the future of 




FIG. 2. A schematic representation of the spacetime struc- 
ture of a realistic black hole formed by the collapse of a ro- 
tating star. An observer falling in with the surface of the star 
encounters a spacelike singularity indicated by the zig-zag line 
in the figure. Gravitational radiation, originating from the 
stellar collapse, gets partially scattered into the black hole; 
at late times, this flux decays as an inverse power of external 
advanced time v~ q . The radiation which crosses the event 
horizon (EH) gets scattered once again by the inner gravita- 
tional potential barrier (indicated in the figure). The result is 
a flux of gravitational waves irradiating the Cauchy horizon 
(CH) and causing it to contract. Note that CH is located 
at v = oo in these coordinates. The inset shows the set up 
of the characteristic initial value problem solved in the text. 
The two initial characteristics E + : u = uo and E~: v = vq 
are indicated. The shear along these characteristic surfaces 
has an inverse power-law dependence on the advanced and 
retarded times reflecting the behavior of the late time wave 
tail of gravitational collapse. 

The late-time wave tail of gravitational collapse pro- 
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duces a flux of radiation across the black-hole event hori- 
zon which decays as an inverse power-law of external 
advanced time p9| . Starting from these initial condi- 
tions Ori has examined the scattering of a test field 
inside a Kerr black hole. His results indicate that the 
amplitudes of test fields decay as an inverse power-law 
in both retarded and advanced times u and v near to 
V in Fig. H Ori also argues that the decay of non- 
axisymmetric modes is modulated by oscillatory terms 
which originate from the rotation of the inner horizon 
with respect to infinity, i.e. the oscillations are a di- 
rect consequence of frame dragging in the Kerr geometry. 
Since these oscillations have no counterpart in spherical 
symmetry it seems worthwhile to outline Ori's argument 
here. 

The line element for the Kerr spacetime, written in 
familiar Boyer-Lindquist coordinates, is 

ds 2 = -(1 - 2Mr/p 2 )dt 2 + (p 2 /A)dr 2 + p 2 d9 2 
+ [r 2 + a 2 + (2Mra 2 /p 2 ) sin 2 9} sin 2 Od(j) 2 
-{AMra/p 2 ) sin 2 Odcpdt , (3.1) 

where p 2 = r 2 + a 2 cos 2 and A = r 2 — 2M r + a 2 . The 
solution depends on two parameters: the mass M, and 
the angular momentum per unit mass a. The equations 
governing perturbations of the Kerr spacetime are not 
fully separable in coordinate space^, this indicates that 
the evolution of various multipoles are coupled. However, 
Ori argues that the coupling between multipoles should 
be weak at late times so that a good first approximation 
to the fields can be obtained by decomposing them over 
the spherical harmonics, and solving the resulting equa- 
tions while neglecting the coupling completely. One then 
iterates to obtain better approximations to the solutions. 
At late times, on a surface of constant coordinate r out- 
side the event horizon, the field is given approximately 
by 



F v {6)e imn - v v- {2l+ V + F u {e)e lmn - U u 



(21+2) 



lm 



F(r, e)e im4 't- { - 2l+2 ^ 



(3.2) 



where I > |m|, and F is some function of both r and 
0. The time dependence of this result has been verified 
numerically by Krivan et al |31|. Now, the coordinate <f> 
is badly behaved at the black hole event horizon. When 
the field is expressed in terms of the regular coordinate 
4> + = (j> — ft+t, where fl + = a/(2Mr + ), and matched 
to the ingoing solution at the event horizon r = r+, one 
arrives at 



F + {9)e im( ^ + e imn + v v -( 2l + 2 ) 



(3.3) 



In the final step, this ingoing solution gets matched to 
the solution near to the Cauchy horizon (denoted by r_ ) 
giving the final expression 



(3.4) 



Here F u and F v are functions to be determined, and 
<f>- = 4> ~ fi-t- The source of the oscillations in the 
non-axisymmetric modes is now obvious. 

While this argument seems plausible, it is far from rig- 
orous. In order to assess the complete significance of 
these non-axisymmetric modes it would be necessary to 
establish the amplitudes of these terms at late times, in- 
formation which is currently unavailable. For this reason 
we focus our attention on the power-law decay of the 
wave tail; however, we do indicate where these oscilla- 
tions might modify the analysis. 



B. Coordinate conditions 

While the conclusions of our analysis are couched in 
terms of physical observables, such as the tidal forces 
(and distortion) experienced by observers approaching 
the singularity inside a black hole, or in terms of curva- 
ture scalars, it is extremely important to understand the 
coordinates used to describe the spacetime. 

The Cauchy horizon can be thought of as the extension 
of future null infinity inside the black hole; that is, the 
Cauchy horizon is located at infinite external advanced 
time. Therefore it is convenient to fix the coordinate v 
to be standard, external advanced time as measured by 
an observer far outside the black hole. Since the exter- 
nal geometry settles down to a stationary state at late 
times, and the strength of the tail radiation crossing the 
black hole horizon decays, we assume that e A ~ e~ KV as 
v — > oo; this is known to hold in non- linear evolution of 
scalar fields in the spherical case, and is the coordinate 
expression of the infinite gravitational blueshift between 
the external universe and the Cauchy horizon. In the 
subsequent analysis we will show that this assumption 
leads to a self-consistent picture of the black-hole inte- 
rior. The approach adopted here is similar in spirit to 
that of BKL S ; we neglect terms in the field equations 
which are suppressed by the exponential factor e A and 
solve the resulting asymptotic equations. 

The coordinate u is chosen so that it goes to negative 
infinity at the event horizon of the black hole. In this way, 
it can be taken to coincide with the natural retarded time 
coordinate inside a Kerr black hole at late times. 

Finally, the coordinates 6 a on the two surfaces of foli- 
ation will be fixed as required by the analysis. 



IV. A SIMPLIFIED MODEL— ALMOST PLANE 
SYMMETRY 



^However, one can separate the equations by Fourier trans- 
forming the fields with respect to t and then completing the 
separation. 



The goal of the present work is to examine the struc- 
ture of spacetime in the neighborhood of the Cauchy 
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horizon of a black hole formed by the collapse of a rotat- 
ing star. The essential features of our analysis are most 
clearly illustrated in a slightly simplified context where 
the Cauchy horizon is irradiated by gravitational waves 
of a single polarization — this interpretation is motivated 
by the local "almost plane symmetry" of the spacetime 
we consider below. 

In this section we set 5^4 = 0, writing the line element 

as 



ds 2 = -2e x dudv + p 2 h ah dB a dB h 



(4.1) 



where we allow the remaining metric functions p 2 , A and 
h a b to depend on all the coordinates {u 7 v,9 a }. In ad- 
dition, we assume that one may simultaneously choose 
coordinates in which the conformal metric h a b is diago- 
nal, and write 



h ab d6 a d6 b = e~ 20 dx 2 



e^dy 2 



(4.2) 



The assumption of vanishing shift vectors reduces the Lie 
derivative operators Da to simple partial derivatives, i.e. 
Da — 9a- As a result, the shear tensor has the simple 
form 



OAab = p 2 /? ^diag[- 



(4.3) 



where a comma indicates partial differentiation. Simi- 
larly, the expansion rates of the null rays orthogonal to 
the surfaces S are simply Ka = 2p~ 1 p,A, and the twist, 



given by Eq. (2.13), vanishes identically. 



Th ese additional assumptions reduce the Eqs. ( 2.14 )- 
(2.16) to a manageable form in which the central fea- 
tures of our arguments are easily understood. Moreover, 
relaxing these conditions in Sec. [v| results in equations 
which are sufficiently similar in structure that only slight 
modification of the following arguments are required to 
complete the general analysis. 



vL,b - \d a K A - -d a D A X + -K A d a X =s , (4.6) 

where we have discarded terms in the equations which 
are damped by a pre- factor e A . 



In the following argument Eq. (4.6) constrains the de- 
pendence of each of the dynamical variables on 6 a near to 
the Cauchy h oriz on; it can be ignored until later. Con- 
tracting Eq. (4.5) with h ab and substituting for Ka and 
Da gives 



(P 2 ),uv - 

which is readily solved for p 2 







P 



+ L(v,e a ) + R{u,e a ) 



(4.7) 



(4.8) 



Here po(0 a ) is assumed to be non- vanishing almost ev- 
erywhere. The functions L(v,9 a ) and R(u,9 a ) are de- 
termined by the initial data. The precise form depends 
on the detailed evolution of the gravitational field be- 
tween the event horizon and E; we make some minimal 
assumptions below. 
Substituting Eq. 



7) into Eq. 



provides a sim- 
ple equation for the shear, which can be expanded to a 
second order equation for as 



pfl P,v + P,vP,u — ■ 

Defining a function A by 

A = A + In p , 
and contracting Eq. (fO) with tjab gives 



A 



(4.9) 



(4.10) 



(4.11) 



where the inhomogeneous term (J ua b&v ab — 2/3 u /3„ for 
the metric in Eq. (4.2) 



A. The asymptotic form of the equations 



B. Initial conditions 



In the coordinates described above, the asymptotic 
regime of interest is characterized by the exponentially 
small value of e A . This is so because we have tai- 
lored our coordinates to the physical mechanism which 
underlies the Cauchy horizon instability — the gravita- 
tional bluc sh ift. Moreover, the vacuum Einstein equa- 
tions fl2.14D - (|2 16| ) reduce to 



D (A K B ) o 



1 



+K {A D B} \ - —tjab {D* + K») D E \ ~ 
-±p 2 h ab (D A K A + K A K A ) 

-D A a A ab + ^A(a a b)d - 



(4.4) 
(4.5) 



and 



Initial data for Eqs. (4.4)-(4.6) is provided on a pair of 
characteristic surfaces, one of which crosses the Cauchy 
horizon as shown in Fig. ^|. This is not completely satis- 
factory, as it assumes the existence of the Cauchy hori- 
zon, at least on one initial characteristic surface. Ideally 
we would like to place initial data on a spacelike surface, 
such as E, below the inner potential barrier, but on the 
past boundary of the region of high blueshift. Unfortu- 
nately, this program is impractical in the present con- 
text. Instead we insure that the solution we construct 
is consistent with such initial data — the solution could 
be smoothly matched to that of a stationary black hole 
background with small perturbations. 

The choice of external advanced time v to describe 
spacetime near the Cauchy horizon implies that e A — > 
on the Cauchy horizon. In particular, on the initial 
characteristic surface E+ (u = u ) we write 
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A|+~ A (9 a ) - k(v - v ) + 



(4.12) 



C. The evolution 



where Ao(9 a ) is its value on the two surface S° — 
{E + HE - }, k is a constant with the dimension of inverse 
length, and terms which vanish in the limit of v — > oo 
are indicated by dots. A similar condition is demanded 
on £~ (v = vq): 



A|_ ~ A (9 a ) - k(u - u ) + 



(4.13) 



Up to the surface £ on which effects of gravitational 
blueshifting begin to take hold, the evolution of the de- 
caying wave tail of gravitational collapse is well described 
by perturbation theory. To the future of this surface the 
large blueshift suggests that geometric optics is valid, 
and that ingoing gravitational waves will be negligibly 
scattered by the gravitational field. Consistent with this 
picture we write the initial data for the shear as 



(4.14) 
(4.15) 



where p,(9 a ) and v(9 a ) are unspecified functions on the 
two surfaces which foliate the two initial characteristic 
hypersurfaces. The inverse power-law decay is motivated 
by the behavior of the tail radiation crossing the event 
horizon of the black hole. The precise nature of the out- 
flux is irrelevant, although it must decay sufficiently fast 
as u — > — oo [0. While perturbation theory suggests 
that p = q, we allow for more general behavior. The qual- 
itative picture which emerges below remains unchanged 
provided decays at least as fast as t>~ 3 / 2 but not 

faster than e~ KV . In particular, oscillatory terms which 
modulate the power-law decay, due to differential rota- 
tion of the event horizon and the Cauchy horizon of a 
rotating black hole, do not change the important fea- 
tures of the analysis below. These oscillatory term s wil l 
show up as the higher order corrections in Eqs. ( 4.14 ) 
and ( p~15D . 

Using Eqs. (4.14) and ( 4.15| ) we can now solve Ray- 
chauduri's equation for p 2 on S + (obtained from Eq. (4.4) 
by setting B = A = v) 



(p 2 ), vv -(p 2 ), v d v A = -2p 2 (/3, v ) 2 



(4.16) 



This determines the function L(v, 9 a ), which can be writ- 
ten as 

L(v, 9") * (^)- g+l 2 ^ ( rjf I 1 + °(V«)] (4-17) 

in the large v limit. The explicit form of R(u, 9 a ) is 
unnecessary. It is sufficient to note that d u R(u,9 a ) — ► 
as u — > — oo. 

With this information, it is straightforward to check 
that the curvature diverges as v — > oo along E + , and that 
the behavior is consistent with the discussion in Brady 
and Chambers Ell]. 



We now examine the solution determined by the initial 
data constructed in the previous section. The evolution 
of the gravitational degrees of freedom (the conforms! 
two-metric h a b) is dictated by the linear wave Eq. (p~S|). 
Therefore the crucial step in determining the geometry 
of spacetime near the Cauchy horizon is to understand 
the evolution of (3 from which the gravitational shear is 
directly determined. If [3 was to diverge the assumptions 
stated in Sec. [II would be violated, and our entire anal- 



ysis would break down. 

In this almost plane symmetric m odel, we can con- 
struct an explicit solution of Eq. (4.9) and demonstrate 
that j3 remains finite, and use it to show that e~ A — > 
on the Cauchy horizon for a finite range of u. For the 
general case presented in Sec. |v| we are not afforded the 
luxury of an exact solution, therefore we also present a 
method which provides a useful bound on the shear near 
to the Cauchy horizon singularity and is applicable in the 
general case. 



The linearity of Eq. (4.9) makes it directly amenable to 
Fourier techniques (see Appendix , however it is more 
illuminating to write the solution as a series 

p = O + p- 1 A lP - 2% [><-*> {V) + G { -^ (U)] (4.18) 



where V = L(v,9 a ), and U = R(u,9 a ). The nota- 
tion F(~^ indicates the i-th integral of th e fun ction 
F(x) = F(°\x) with respect to x. Inserting (4.18) into 



Eq. (4.15) and comparing integrals of the same order 
gives the coefficients 



At 



A 



(2z + l) 2 
4(i + l) 



(4.19) 



The functions F(V), G(U) and A a (9 a ) are determined 
by the initial data. 



Differentiating Eq. ( 4.1Sj ) with respect to v and com- 
paring with the initial data along S + determines F in 
the limit v — > oo to be given by 



ZAq 



(4.20) 



from which one can extract the behavior of F to be 



F(v, 



A k(2 - q) 



(4.21) 



Given the solution ( |4.2l[ ) and L(v, 9 a ) in Eq. ( J4.17 ) one 
derives the recursion relation 



pi-i-i) — p(-i). 



(kv) 
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— . (4.22) 
(i + l)(l-«) + l-g/2 1 ; 

Combining Eqs. fl4.22| ), ( f4.21| ) and ( [1.19[ ) it is easy to 
check that the sum over converges to a finite re- 

sult provided q > 3 — the boundary conditions for the 
radiative tail at the event horizon actually imply q > 14. 







It follows that the leading order term in Eq. ( |4.1£j ) is 
proportional to F since F ( -^(V) > F < -- l -^(V) for large 
v. This shows that /3 „ continues to decay with an inverse 
power-law form all along the Cauchy horizon provided it 
has such a decay on the initial surface. Similar arguments 
provide G. 



We c an no w use Eq. (4. IS) to compute the source term 
in Eq. ( 4.11 ) and hence determine the evolution of A to 
be 



A~ A\-(u,O a ) - k{v - v ) - / du' dv'a uah a v 

Juq Jv 



a b 



(4.23) 



(KV 



-9/2 



The integrand in Eq. (4.23) is proportional to /3,u/3. 



as v — > oo, consequently the integral is finite in 
this limit provided q > 2 as we have already required. 
Clearly, A diverges to negative infinity at the Cauchy 
horizon in precisely the same manner as it does along 
the initial hypersurface X + , that is 



-KV 



(4.24) 



) provides a final consistency 



A 



as v — ► oo. Equation 
check on the solution. 

Ha ving dete rmined the approximate solution to 
Eqs. (4.4)-(4.6) near to the Cauchy horizon we are in 



a position to examine the curvature. It is convenient to 
consider the Newman-Penrose components of the Weyl 
tensor on the basis 

{e- A/2 %,), e- x 'H (1)l e (o) m°, e (a) m a } , (4.25) 

where m a is a complex two- vector (called the shear axis) 
which satisfies g a b = 2m( a ™>b)- The leading behavior of 
each of the Weyl scalars is presented below: 

^>u=e- x (Kv)-i[-KLi{e a ) + ...} 



*2= e- x (nv)-i [-ii{O a )v(O a ){-Ku)-% + . . 
*3= (nv)^^ +1 [fh a d a (iJ,v)(— ku)~% + . . 
* 4 = e- x (-Ku)-% [-Kv{8 a ) + ...}. 



(4.26) 



Notice that M^i and ^3 are non-zero, this arises because 
the spacetime is not exactly plane symmetric. The square 
of the Weyl curvature C a /3 1 sC al3l ' s diverges as v — > 00 at 
the Cauchy horizon, and is dominated by the radiative 
piece ^0^4 on this tetrad. Degeneracy of the roots of 
the polynomial 



a * 4 + ar<3> 3 + a * 2 + a*i + *o = 



(4.27) 



determines the Petrov classification of the spacetime [g8| . 
Brady and Chambers [|T) demonstrated a four-fold de- 
generacy in the limit that v — > 00 on S + . Using 
Eqs. (4.26) one shows that this statement continues to 
hold near to the Cauchy horizon but away from E + , that 



is, the diverging curvature is asymptotically Petrov type 
N. This suggests the intuitive physical picture of the sin- 
gularity as a gravitational shock wave propagating along 
the Cauchy horizonl. It is worth noting however that 
this approach to type N behavior is not characterized by 
a peeling property as it is at large distances outside the 
black hole. Indeed a 2 — > ±\/— ^0/^4 as f — > 00 so that 
all four roots tend to zero. This behavior is indicated 
schematically in Fig |[ 



Im 





Re 



FIG. 3. The complex plane showing a schematic represen- 
tation of the roots of Eq. ( 4.27 ). The arrows indicate the evo- 
lution of the roots as v — * 00. The roots become degenerate 
in the asymptotic limit, but they do not exhibit any peeling 
like properties. This suggests that the Cauchy horizon singu- 
larity is Petrov Type N, and can be thought of intuitively as 
a shock wave propagating into the black hole. 



D. Alternative bounds 

The most important step in validating the approxi- 
mation adopted in this paper is to show that e A — > 
as v — > 00 everywhere along the Cauchy horizon. In the 
simple model s pacet ime examined above, we have the ex- 
plicit solution ( 4.18 ) for the metric func tion (3 which de- 
termines the shear tensors through (4.3). Once the shear 
tensors arc known, the evolution of A is given by the inte- 
gral equation (4.23). The linearity of the wave equation 
(4.9) allows the direct computation of/?, and t he ex plicit 
verification that the integral appearing in Eq. (|4.23| ) van- 
ishes in the limit v — > 00. 

In a general spacetime, the evolution of the shear ten- 
sors are described by a non-linear wave equation for 
which no exact solution is known. However, it is still 
possible to place bounds on the integral appearing in 



*This interpretation is somewhat simplistic since the curva- 
ture scalar Cctf3-,sC al3 ~' s vanishes identically in the case of a 
gravitational shock wave. 
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Eq. ( 4 . 23| ) by using the field equations, thus verifying 
that e 7 ^ as v — > oo. Although this procedure is un- 
necessary here, it is instructive to apply it first to this 
case before considering the general non-linear problem in 
the succeeding section. 

First define a new function £ which is proportional to 
the integrand, 



0°) = 2/3 lU /3,, 



(4.28) 



Multiplying the wave Eq (4.9) for 8 by 3 }V , it is not 
difficult to derive 



(4.29) 



Integrating ( 4.29p with res pect to u, making use of the 
characteristic initial data ( |4. 14 ) and the solution (4.17) 
for the v dependence of p 2 , we find 



np 2 (u, v, 9 a ) 



n ; / ) (4.30) 

"0 



for u > u . 

Suppose that at a point with coordinates u, v where 
u € [uq, u] and v £ [vq,v], the largest value of the function 
£ occurs such that £(tt, v) = £. We denote the value of 
any function evaluated at this point with an overbar. The 
partial derivative of 3 evaluated at this point is bounded 
by 



[p{e a )p {6 a )Y(Kv)-i 

up 2 



(/c+(u-tio)e). (4.31) 



A similar bound on the square of the partial derivative 
of 3 with respect to u can be derived, 



3 2 < 



nnpo(n?(~™)- p 

up 2 



(«+(tJ-«0)e). (4.32) 



It then follows from these bounds and from the definition 
of £ that the inequality 

f - 2 J% 

( Pvpl V 



: 4 l^^rj {Kvy q {-nu)-v (k + (u - u )£) 



(4.33) 



must be satisfied if the field equations are satisfied. This 
inequality can be rearranged into the form 



where the coefficients are given by 



(4.34) 



a = l 



\ kp 2 J 



(Kv)- q (-Ku)- p {u-u )(v-v ) (4.35) 



b = 4 K (^4-] ( K vy q (-Ku)- p (u-u o +v-v o )>0 (4.36) 
V K P J 



c = 4k 2 (Wo) ( K v)-i(- K u)-p>0. 
V K P J 



(4.37) 



Clearly, if the coefficient a in this inequality is negative, 
then no bound can be placed on the maximum value of 
£. However, if a is positive, then the upper bound on £ 
is just given by £ < b/2a + (b 2 + Aac) 1 / 2 /2a. Hence, it 
is important to determine the sign of the coefficient a, 
which will depend on the relative magnitudes of the two 
contributions to a. 

Now, consider the magnitude of the second term in 
(4.35). At a fixed value of u, this term is proportional 
to (nv)~ q n(v — vq). In the characteristic diamond (the 
region above the characteristics E + and S~ in Fig. ||) the 
advanced time coordinates v and vo obey 

KV > KVq ^> 1. 

As a result, (nv)~ q n(v — vq) <C 1 in the region of interest. 
The u dependence of this term is \ku\~ p n(u — Uq). Note 
that for all points in the characteristic diamond, < v 
and | it | < |«o|- It then follows that 



K,(u — Uq) < k\uq\ < KV 

so the second term of a is proportional to 

k(u — uq)k(v — vo){nv)~ q \Ku\~ p < (nv)~~ q+2 \Ku\ 



- p « 1, 
(4.38) 



if q > 3. The prefactor 4 ^ ^ which multiplies the ex- 
pression Q4.38[ ) in Eq. ( |4.35| ) is finite, and as a result, the 
coefficient a is positive and approximately unity. Hence 
the quadratic inequality can be used to place the follow- 
ing limit on £, 



e<2^t M 
P 



-q/2 



(—KU) 



-p/2 



(4.39) 



A recent calculation j25) suggests that the corre ct ini - 
tial data for 3 should be an inverse power law (4.14) 
modulated by an oscillatory function of the form cos(kw). 
Note that the arguments leading to the bound ( 4.39| ) will 
be unaltered by a modulation of this sort, since the cosine 
function is bounded by o ne. 

Given the bound (4.39) on £, we are now in a position 
to integrate Eq. ( 4.23 ) and solve for the metric function 
A, 

A(u,v,d a ) < A|_(u,6» a ) - k(v-v q ) - (u-u )(v - v )(. 

(4.40) 

Since we have shown that £ is bounde d by a vanishingly 
small function, the integrated term in ( 4.40 ) is negligible 
compared to the homogeneous terms which arise from the 
initial data. Hence, e A vanishes at the Cauchy horizon 
for all u > Uq. 
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V. THE GENERIC CASE 

Our attention so far has focused on the toy model of 
the interior presented in the previous section. Our ap- 
proach was to use a metric with only three free functions 
of four variables, and choose initial data motivated by 
the theory of scattered fields on a stationary background. 
Although the almost plane symmetric model of the previ- 
ous section does not have sufficient degrees of freedom to 
describe the evolution of generic gravitational perturba- 
tions, the model is useful since we were able to write an 
explicit solution and show that a weak, null curvature sin- 
gularity occurs at the spacetime's Cauchy horizon. The 
importance of the toy model will become apparent in this 
section where we show that the metric of a general space- 
time near to a Cauchy horizon is, to leading order, nearly 
identical to the almost plane symmetric metric. 

The metric used to study the generic evolution of grav- 
itational perturbations is the general 2 + 2 metric (2.9). 
As this metric has eight functions, there exists the gauge 
freedom to set two functions to zero. We choose to set 
the shift vector 8® = 0. As a result, the normal Lie 
derivative operator D v reduces to the partial derivative 
d v . The general metric which we use is then, 

ds 2 = -2e A du dv + p 2 h ab {dd a + s a u du){dO b + s b u du). 

(5.1) 

where the conformal two-metric h a b has unit determinant 
h = det||/i a f,|| = 1 and hence only represents two free 
functions. In the previous section we set one of these 
functions to zero, but here we will consider the evolution 
of the general form of the conformal two-metric. For the 
present, we will not choose any particular representation 
for h a b- 

The choice of coordinates and characteristic initial 
data are motivated by the discussion presented in Sec. |TJ. 
As a result of our coordinate conditions, we assume that 
on the initial characteristics E + and £~ the f orm o f the 
function A = A + In p is given by equations ( 4. 12| ) and 



( 4.13 ) as in the plane symmetric spacetime. The power 
law initial data for the gravitational perturbations can 
be set by specifying h a b on the initial characteristics or 
equivalently, by setting the values of the shear tensors to 



KV 



0~vab\ + — V2( 

a uab \- ~ v2(- 



-9/2 



MO") H 
Vab{B a ) 



-p/2 



(5.2) 
(5.3) 



where p, a b and are traceless two-tensors, h ab u a i, = 
h ab v a b = and the shear tensors are defined by GAab — 
\p 2 DAh a b- As in the previous section oscillatory terms 
have been neglected in the initial data because their pres- 
ence does not alter the qualitative features of the subse- 
quent analysis. 

The analysis of the generic case is performed in the 
same spirit as in the "almost plane" case, and follows 
closely the steps taken in the last section. As before, the 



solution of the metric function p on a u = const, charac- 
teristic is found by solving Raychaudhuri's equation, 

(P 2 ),w - (p 2 ),v cU = VK| 2 , (5.4) 

where we have defined the norm of a two-tensor to be 

\&a\ 2 = OAaOAb" ■ (5.5) 



Given the assumptions ( 4.12| ) and ( |5.2[ ) for the behaviors 
of A and a va b on S + , we find asymptotically, 



p 2 (9 a ) + L(v,9 a ) 



(5.6) 



where L{v,9 a ) is the same function ( 4.17 ) found in the 
plane symmetric spacetime. 



A. The shift and twist vectors 

The non-vanishing shift means that the congruence of 
null rays undergoes some twist. Given our gauge choice, 
the twist is related to the shift by 



uj a = -d v s a u 



(5.7) 



The vacuum field equation R va = ( p.lq ) specifies the 
behavior of the twist on a constant u hypersurface: 



A -a 



{d v +K V -d v X)u a = e A j 



(5.8) 



where 



f = 2(a * 6 ~ kd a K v ~ \d a d v \ + \K v d a \) (5.9) 



is independent of u> a . We can formally integrate Eq. ( |5.8| ) 
to get 



(f 



u a (u,9 a ) 



dv'p 2 f 



(5.10) 



where the function u>q(u, 9 a ) is determined by the initial 
data on S". On S + , the v dependence of the functions 
j a \ + arc known to be inverse power laws, so the integral 
in ( 5.10| ) vanishes asymptotically and the twist behaves 
like 

e A 

u a \ + = -^uj$(u , 9 a ) -> exp{- K v) , v -> oo . (5.11) 
P 

Hence the twist vector is exponentially suppressed on S + . 

In our problem, we still have the freedom to fix co- 
ordinates along one hypersurface of constant v p4|] . A 
natural choice is to ask that given coordinates 9 a at one 
point on the Cauchy horizon, they stay constant when 
Lie convected by £ u along the Cauchy horizon. This is 
equivalent to the statement that s" = on the Cauchy 
horizon. Given this boundary condition, and a solution 
for t he tw ist, the shift vector can be found by integrating 
Eq. (5.7). On S + the shift vector behaves as 
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1 4- ~ exp(— kv) , v 



(5.12) with positive definite entries, an d the norm of a two ten 



On the initial characteristic S + the shift vector is also 
exponentially suppressed. 

The behavior of the shift and twist vectors on later 
hypersurfaces will depend on the behavior of the func- 
tions j°, and hence of the shear and the metric functions 
p and A on later hypersurfaces. In the next part of our 
analysis, we will assume that the shift and twist vectors 
are exponentially suppressed on later hypersurfaces, solve 
the resulting equations, and show that this assumption 
is self-consistent. The effect of this assumption is that 
to leading order in e A , the shift and twist vectors drop 
out of the vacuum evolution equations. The field equa- 
tion ^R a bg ab = reduces to equation ( j4.7| ) for p. As 
in the plane symmetric spacetime, the general solution 
for p is p 2 = pi + L(v, 6 a )+R(u, a ). The field equation 



for A on l ater h ypersurfaces is again given by the integral 
equation (1.23). 



B. The evolution of the shear 

The most significant difference between the general 
light-like geometry and the model in section IV is the 
form of the shear tensors, (JAab- In section IV we as- 
sumed that the shear tensors have a simple diag onal form, 
which resulted in a propagation equation (4^) linear in 
shear. As a result, it was fairly simple to show that 
regular power law initial data for the shear evolves via 
the field equations to a regular solution on later hyper- 
surfaces which also decays as an inverse power law. In 
this section we allow a general (non-diagonal) form for 
the shear tensors. Although the resulting wave equation 
describing the evolution of the shear is non linear the ar- 
gument presented in the last section can still be used to 
place a limit on possible divergences of the shear. 

The evolution of the shear tensors is governed by 



1 



= WR\- ^-S a h ^R d 



e- x (D A a A a h + K A a 



A b\ 



0(1), 



(5.13) 



where we have made use of the fact that the twist and 
shift vectors are exponentially sup press ed in our approx- 
imation scheme. Although Eq. ( 5.13 ) is linear in the 
shear tensor, it should be noted that <JAa = \h ac dAh ac 
which is n on- linear in the conformal two-metric. Thus 
Eq. ( 5.13| ) is a non- linear wave equation for the two- 
metric. 

It is useful to introduce the following matrix 



, 1 / & a l b D a 



= g- dia s(l cr «l 2 >KI 



(5.14) 



sor 
is 



\<Ja\ was defined in Eq. ( p.5|) . The divergence of tAB 



d BtA B = — {(TBb a d A (T A a b + 2a A b a d [A a B]a b ) . (5.15) 

The anti-symmetric derivative of the shear occurring in 
the second term of ( 5.1 5| ) is related to the twist as has 
been shown in appendix B of reference BJ] and is thus 



exponentially small compared to th e first term of ( 5T5). 
Substituting the vacuum field Eq ( ^13| ) into (| .151) we 
find the following evolution law for the components of 

tAB'- 

\d B (p 2 t A B ) = —^-K A ° D a b o D b a + 0(e x ). (5.16) 



R uv = reduces to Eq. (|4.1l| ) for A. The formal solution The quantity 



16tt 



C ~ b „L> a „ 

4 = — ^ Da b = 



b a 
ua Vvb 



(5.17) 



on the right hand side of this equation is the same quan- 
tity appearing in Eq (4.11) for t he m etric function A. By 
placing bounds on £, using Eq f[l|), we can show that 
the functional form of the initial data A|+ on S + is a good 
approximation for the form of A on later hypersurfaces. 

Fir st note that the definition ( |5 . 1 7| ) for £ r educ es to 
(4.28) in the case of the plan e wa ve metric (4.1). In 
component form the field Eqs ( |5.16| ) are 



d u { P 2 \o v \ 2 )= -d vP 2 ii 

d v {p 2 \o u \ 2 )= ~d uP 2 ^ 



which are reminiscent of Eq 
isfies the Schwartz inequality. 



1" 

HE 



(5.18) 
(5.19) 

In fact, since £ sat- 



< \<j v 



the ar- 



gument encapsulated in Eqs (4.29 - 4.40| ) hold for the 



general dou ble-null metric given the assumption made in 
section III. 



C. Solution of the initial value problem 

As mentioned before the dynamical degrees of freedom 
are encoded in the shear tensors. Once we have fixed a 
gauge and know the dynamic evolution for the shear we 
can in principle calculate all the metric functions. Initial 
data for the shear is supplied on £ in the form of a pair 
of traceless tensors cr va b(v, 9 a )\ + and a ua b{u,& a )\- given 
by Eqs. (^2) and ( |5.3| ). Since the evolution equations for 
the shear are non-linear we do not expect to find a closed 
form solution, however we can cons truct an approximate 
solution along the lines of Eq. ( 4.18 ). Write the conformal 
metric explicitly as 



, _ / e 2/3 cosh 7 sinh 7 \ 
ab ~ \ sinh 7 e 2 ? cosh 7 ) ' 

The equations for the functions (3 and 7 are then 



(5.20) 
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[p cosh(7)y3 )V ] 



f [pcosh(j)). v (3^ 



(5. 

2sinh(2 7 )# l ,/?„ , (5. 



where we have neglected the terms involving the shift 
which is exponentially suppressed according to the argu- 
ments above. It is reasonable to assume that (p cosh 7)^ 
should be weak ly dependent on u, so that we can inte- 
grate Eq. (5.21) for (3 and obtain 



21) weak in the sense that a coordinate system exists in which 

22) the spacetime metric is regular at the singularity A 
similar result holds in the context of the approximate so- 
lution presented in the previous section. By introducing 
a new coordinate V = —e~ KV , the asymptotic form of the 
line element is 



B hl (v,x a )+B out (u,x a ) 
p cosh 7 



(5.23) 



Notice that this reduces to the leading term in the se- 
ries solution for (3 presented in Eq. (4.18) when 7 = 0. 
Since the shear is bounded and small according to the 



arguments in subsection VB , we further expect that the 
non-linear term in Eq. (5.22) can be treated as a source, 
that is, we assume that 7 is slowly varying in the region 
of interest to us. Moreover, this term involves the prod- 
uct (3^ v (3 iU which is effectively quadratic in the luminosity 
of the infalling gravitational wave tail, and therefore less 
important than the boundary terms. Thus we have 



7 



G in (v,x a ) + G out (u,x a ) 



(5.24) 



Finally, we fix the four free functions which appe ar in 
this solu tion by reference to the initial data in Eqs. (5.5) 
and (5.2). The validity of these approximations has been 
confirmed by Droz by numerically integrating the equa- 
tions Q. 

For vanishing shifts Eqs. ( |4.4| )-(4.6) continue to hold 
as no assumptions about the form of the shear have been 
made in their derivation. The analysis proceeds along 
the same lines as in the previo us se ction, except that the 
source to the wave-equation (4.11) has a more compli- 
cated functional form, but it is still small and decaying. 
We therefore recover the result in Eq. (4.24) for A. Sim- 
ilarly the function p 2 is recovered from integrating Ray- 
chauduri's equations. 

It is now straightforward to check that the asymptotic 
beha vior o f the Weyl scalars is of the same form as in 
Eqs. ( 4.26 ) — the terms inside the square brackets are dif- 
ferent for the general case, however the scaling in u and v 
arc identical. Hence the intuitive picture of the singular- 
ity as a gravitational shock propagating along the Cauchy 
horizon continues to be valid in the generic case. Thus, 
we have demonstrated that the generic structure of the 
Cauchy horizon singularity is qualitatively cap tured by 
the almost plane symmetric model of section IV , and the 
Cauchy horizon singularity occurs for generic perturba- 
tions (as long as the initial data is not too singular |32| ) . 



VI. STRENGTH OF THE SINGULARITY 

One of remarkable things about the mass-inflation sin- 
gularity inside charged, spherical black holes is that it is 



ds 2 



2e A 

KpV 



dudV + p 2 h ab (dd a + s a u du){d6 b + s b u du) 



(6.1) 



Combining Eqs. ( 4.24[ ) and ( |4~g| ) with the definition of V 
it is clear that ( 4 >g a /3 is bounded as V — ► 0, i.e. at the 
Cauchy horizon. 

What is the relevance of this coordinate system to ob- 
servations? This is somewhat clarified by re-stating the 
result; there exists a coordinate system in which twice in- 
tegrating the curvature with respect to the new advanced 
time V gives a finite result. It turns out that the proper 
time t measured by an observer crossing the Cauchy hori- 
zon satisfies r ~ V, so the tidal acceleration experienced 
by this observer diverges like r _2 (ln |r|)~™ as the null 
singularity is approached, where n > 21 + 3. The tidal 
distortion is given by twice integrating this acceleration 
along the worldline of the observer, and it is finite all 
the way up to the singularity. This rough argument sug- 
gests that the regularity of the metric, when written in 
terms of the coordinate V, indicates that the singularity 
at V = is weak, in the sense that tidal distortion of an 
extended object is finite there. 

The interpretation of this result is somewhat unclear. 
One might be tempted to think about the singularity 
along the Cauchy horizon as an "impulsive" singularity — 
while an infinite force is exerted, it acts only for a very 
short time. Such a viewpoint has been adopted by some 
authors and taken to indicate that a classical continua- 
tion of spacetime beyond the Cauchy horizon singularity 
might exist [Q,^J . Unfortunately, this point of view seems 
problematic since classical physics provides no mecha- 
nism by which to regulate the curvature once it diverges: 
only pure gravitational shock waves can be confined to a 
thin layer in classical General Relativity |Q. Indeed, we 
know that quantum effects are important in the descrip- 
tion of spacetime near the Cauchy horizon singularity 
(see, for example, the discussions in 1 33 34]), and may 



dramatically change the classical picture. 



VII. CONCLUSION 

Our results strongly indicate that a wide class of initial 
data can lead to the formation of a weak, null curvature 
singularity inside a black hole formed by gravitational 
collapse. The analysis, which is valid at late times near 
to the singularity, demonstrates that the null character of 
the singularity is independent of the initial data provided 
the flux of radiation entering the black hole at late times 
falls off more quickly that v~ 3 , that is the shear o va bVv ah 



14 



decays at least this fast along the event horizon. More- 
over, our approximate solution depends on 2 x 2 func- 
tions {a ua b(u ,v), cr va b{u,Vo), traAab = 0) corresponding 
to the physical degrees of freedom of the gravitational 
field, and we can therefore claim that the null singularity 
is generic and not an artifact of special symmetry. 

It is important to compare our results to those of previ- 
ous analyses. Ori has investigated the singularity inside 
a realistic, rotating black hole using non-linear pertur- 
bation theory ||. The picture of a weak, null singu- 
larity that emerges from our work is in agreement with 
the results of his analysis. Furthermore, our results lend 
support to Ori's claim that the perturbative approach 
captures the essential features of spacetime structure in 
the neighborhood of the Cauchy horizon singularity. 

Some work remains to be done, however. In spherical 
models, the null Cauchy horizon singularity is a precursor 
of a strong spacelike singularity deep inside the black hole 
core A similar result presumably holds for real- 

istic rotating black holes, however little is kn own a bout 
this situation. We have also seen in section III A that 



linearized perturbations of the Kerr black hole may be 
modulated by terms which oscillate infinitely many times 
as v — > oo; it is important to determine how significant 
these oscillations are for the variation of curvature as 
measured by an observer approaching the Cauchy hori- 
zon singularity. Two approaches seem worth pursuing to 
further explore this issue: (i) linear and non-linear per- 
turbation theory |l^J^| can provide an answer, in princi- 
ple, however it requires the difficult computation of the 
relative amplitudes of all the terms in the perturbation 
series; (ii) An alternative approach is provided by nu- 
merical techniques, similar to those used in the spherical 
case Ip^f . While double null formulations of numerical 
relativity encounter serious problems when caustics form 
along the characteristic surfaces used in the evolution, we 
have in this problem a well understood regime in which 
such coordinate difficulties can surely be overcome. 
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APPENDIX A: EXACT PLANE-SYMMETRIC 
SOLUTION 

Yurtsever has proved that singularities which form in 
colliding plane wave spacetimes are generically space- 
like |35j. In this context, the exact solution to Eq. ( 4~9| ) 
was also worked out by Yurtsever p5| . For complete- 
ness, we present his derivation here and clarify why the 
isomorphism between the internal geometries of rotating 
black holes and colliding plane wave spacetimes docs not 
imply that singularities inside black holes are generically 
spacelike. See Ref. Q for a related argument. 

Eq. (4.7) for the conformal factor p 2 can be viewed as 
an integrability condition for a coordinate transformation 
from the null coordinates (u, v) to new coordinates (a, x) 
defined by 



a = p 2 = p 2 + L(v,d a ) + R(u,6 a ) 
X = L(v,9 a )-R(u,8 a ) . 



(Al) 
(A2) 



Notice that a is a time coordinate since we are inside 
the event horizon of the black hole. In terms of these 



coordinates the wave equation (O) for (3 becomes 



ft,aa AxX 



a 







(A3) 



Since Eq. (A3) is manifestly independent of time, we 
eliminate x using a Fourier transform. The result is 
Bcsscl's equation in a so that the solution is 

(3 = p + [dke lkx [c(k)Y (\k\aj + d(k)J (\k\a)} (A4) 



where Jo and Yq are Bessel functions of the first and 
second kinds, and the two functions c(k) and d(k) are 
determined by the initial conditions on E in Fig. |2|. The 
function Y"o(|^l a ) diverges as a — ► 0, but since a is always 
larger then zero inside a black hole (at least up to early 
portion of the Cauchy horizon) (3 should remain regular 
for all physically relevant values of a and x inside black 
holes. 

In contrast to the black hole case, the Cauchy horizon 
in colliding plane wave spacetimes is at a = and conse- 
quently the gravitational shear generically diverges there; 
this causes the catastrophic focusing of ingoing lightrays 
to a spacelike singularity. Arguments demonstrating that 
spacelike singularities are generic in colliding plane wave 
spacetimes cannot be directly generalized to the black 
hole interiors. 



[1] S. W. Hawking and G. F. R. Ellis, Large scale structure 
of spacetime (Cambridge University Press, Cambridge, 
1973). 



15 



[2] V. A. Belinsky, I. M. Khalatnikov, and E. M. Lifshitz, 



qc/9801086| ). 



[37] 



Sov. Phys. JETP 32, 169 (1970). 

E. Poisson and W. Israel, Phys. Rev. D 41, 1796 (1990). 

A. Ori, Phys. Rev. Lett. 67, 789 (1991). 

A. Ori, Phys. Rev. Lett. 68, 2117 (1992). 

R. Penrose, in Battelle Rencontres, edited by C. de Witt 

and J. Wheeler (W.A. Benjamin, New York, 1968), p. 

222. 

S. Chandrasekhar and J. Hartle, Proc.R.Soc.Lond. 
A384, 301 (1982). 

R. Price, Phys. Rev. D 5, 2419 (1972). 

C. Gundlach, R. Price, and J. Pullin, Phys. Rev. D 49, 

883 (1994). 

W. Hiscock, Phys. Lett. A 83, 110 (1981). 

R. Isaacson, Phys. Rev. 160, 1272 (1968). 

Y. Gursel, I. D. Novikov, V. D. Sandberg, and A. A. 

Starobinsky, Phys. Rev. D 20, 1260 (1979). 

A. Ori, Phys. Rev. D 55, 4860 (1997). 

A. Bonanno, S. Droz, W. Israel, and S. M. Morsink, Phys. 

Rev. D50, 7372 (1994). 

R. Balbinot, P. R. Brady, W. Israel, and E. Poisson, 
Phys. Lett. B 161, 223 (1991). 

R. Hermann and W. G. Hiscock, Phys. Rev. D 46, 1863 

(1992) . 

M. L. Gnedin and N. Y. Gnedin, Class. Quant. Grav. 10, 
1083 (1993). 

P. R. Brady and J. D. Smith, Phys. Rev. Lett. 75, 1256 
(1995); L. Burko, Phys. Rev. Lett. 79, 4958 (1997). 
L. Burko and A. Ori, "Analytic study of the null singu- 
l arity inside sph erical charged black holes", unpublished 
flgr-qc/9711032| ). 

A. Bonanno, Phys. Rev. D 53, 7373 (1996). 

P. R. Brady and C. M. Chambers, Phys. Rev. D 51, 4177 

(1995). 

O. Ori and E. Flanagan, Phys. Rev. D 53, R1755 (1996). 
U. Yurtsever, Class. Quant. Grav 10, L17 (1993). 
P. R. Brady, S. Droz, W. Israel, and S. M. Morsink, Class. 
Quant. Grav. 13, 2211 (1996). 
A. Ori, Gen. Relativ. Grav. 7, 881 (1997). 
C. Barrabes and W. Israel, Phys. Rev. D 43, 1129 (1991). 
C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravi- 
tation (Freeman, New York, 1973). 

S. Chandrasekhar, The mathematical theory of black 

holes (Oxford University Press, New York, 1983). 

W. Krivan, P. Laguna, and P. Papadopoulos, Phys. Rev. 

D54, 4728 (1996). 

W. Israel, unpublished . 

W. Krivan, P. Laguna, P. Papadopoulos, and N. Ander- 

sson, Phys. Rev. D56, 3395 (1997). 

S. Droz, Phys. Rev. D55, 3575 (1997). 

W. G. Anderson, P. R. Brady, W. Israel, and S. M. 

Morsink, Phys. Rev. Lett. 70, 1041 (1993). 

R. Balbinot and E. Poisson, Phys. Rev. Lett. 70, 13 

(1993) . 

U. Yurtsever, Phys. Rev. D 38, 1706 (1988). 
L. M. Burko, in Internal structure of black holes and 
spacetime singularities, edited by L. M. Burko and A. 
Ori (IOP Publishing, 1998). 

A. Ori, "Null weak singularities in plane-symmetric 
spacetimes", Phys. Rev. D 57, to appear (jr- 



1G 



